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1 Introduction 

The Standard Model (SM) has been compared to experiments with great success in the 
past decades and finding the Higgs boson is the only missing piece. However, there are still 
a few internal problems. The prime example is the so called hierarchy problem: why is 
the electroweak (EW) scale much smaller than the Plank scale? Thus, the SM cannot be 
seen as a fundamental theory of particle physics, but only as an effective description which 
will break down at higher energies, at least at the Planck scale where gravity becomes of 
the same magnitude as the gauge forces. The mission of the Large Hadron Collider at 
CERN is therefore not only to look for the SM Higgs boson but also for physics Beyond 
the Standard Model (BSM). 

The general two Higgs Doublet Model (2HDM) was one of the earliest BSM models, 
proposed by T.D. Lee [1] already in 1973 as a model with spontaneous CP-violation. The 
2HDM itself cannot give any solution to the problems of the SM, such as the hierarchy 
problem. On the contrary, it introduces more problems such as tree level flavour-changing- 
neutral-currents (FCNC) which are absent in the SM. However, a 2HDM is part of many 
other BSM models, especially supersymmetric ones, which require an even number of Higgs 
doublets. Therefore it is useful and interesting to study the 2HDM itself, since it can be 
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thought of as an effective description of more general models at the TeV scale. One such 
example is the Minimal SuperSymmetric Model (MSSM) in the case of heavy superpartners 
such that the Higgs bosons only decays to SM particles. 

The problem of tree level FCNC can be evaded by introducing an appropriate Z2 
symmetry that ensures that each fermion type only couples to one of the Higgs doublets, 
which is sufficient in order to avoid tree-level FCNC as shown by Glashow and Weinberg [2]. 
This is precisely what happens in the MSSM whose Higgs sector at tree-level is a so called 
type II 2HDM, meaning that one of the Higgs doublets couples only to down-type fermions 
and the other only to up- type ones. By enforcing a ^-symmetry one also ensures the 
absence of tree-level FCNC under renormalization group evolution of the model to other 
energy scales. 

Recently another way of avoiding the tree- level FCNC, by having the Yukawa couplings 
to the two Higgs doublets proportional to each other, has been proposed [3]. This works 
fine at a given energy scale but if one evolves the model to another scale then the tree-level 
FCNC are reintroduced because the Yukawa couplings in this model do not respect any 
Z<i symmetry as shown by Ferreira et al [4]. There has also been some discussion of the 
experimental constraints on this model under renormalization group evolution [5-7] and 
we will revisit these constraints more carefully below. 

More generally, the FCNC at a given energy scale are avoided as long as the Yukawa 
couplings are diagonal in the appropriate basis. The constraints on these more general 
models from low-energy flavour observables have also been studied [8], but not their prop- 
erties under renormalization group evolution. Apart from these schemes, which are set 
up in order to avoid tree- level FCNC to a larger or lesser extent, one can also envision a 
top-down approach where one assumes a certain texture for the mass matrices and from 
this derives the Yukawa coupling matrices. In the present context the prime example is 
the Cheng-Sher ansatz [9] which gives a natural suppression of tree-level FCNC from the 
hierarchy of quark masses. Some generic properties of these models under renormaliza- 
tion group evolution have been studied [10] but not taking experimental constraints into 
account. 

In this paper we will study the properties of all these types of models taking into 
account also experimental constraints on FCNC when evolving them according to the 
Renormalization Group Equations (RGE) for the Yukawa couplings. In this way we can 
see how stable the various assumptions are under RGE evolution, which in turn gives 
a measure of have plausible the assumptions are. A large sensitivity indicates that the 
assumptions behind the model are not stable meaning that they are either fine-tuned or 
incomplete such that there for example will be additional particles appearing when going to 
a higher energy. From this respect we will study both the appearance of a Landau pole as 
well as off-diagonal Yukawa couplings leading to FCNC. Strictly speaking, the experimental 
constraints on the latter are given at the EW scale. Even so, we can still apply them at a 
higher scale as a means of a determining the fine-tuning of the models as argued above. As 
an alternative one can also envision to assume a ^-symmetric starting point at the EW 
scale, then evolve up to a high scale where the ^-symmetry is broken, and finally evolve 
down to the EW scale again where the experimental FCNC constraints are then applied. 
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The layout of the paper is as follows. We first give a brief introduction to the general 
2HDM in section 2 including the Yukawa sector with emphasis on the FCNC problem as 
well as some possible solutions and the RGEs for the Yukawa couplings. Section 3 gives 
the latest constraints on the non-diagonal Yukawa couplings from neutral meson mixing as 
well as the SM input values we use. Then in section 4 we present our numerical analysis 
of the running Yukawa couplings. We investigate the limits both from the absence of a 
Landau pole as well as from requiring the off-diagonal Yukawa couplings at higher energy 
scales to be in accordance with the experimental limits at the EW scale. Finally, in section 
5 we present our conclusions. 



2 The general 2HDM 
2.1 The Scalar Sector 

The two Higgs doublet model was introduced in [1] and for a more general overview of its 
properties and the constraints that can be put on it, we refer to the recent review [11]. 
Much of the phenomenology of the 2HDM is also closely related to the SM and MSSM for 
which we refer to the reviews by Djouadi [12, 13]. 

The most general renormalizable scalar potential with two Higgs doublets, $i and 
can be written as 

V$ = mu$}$i + m^^^a - (m? 2 $i*2 + h.c) 

+ 1a 1 ($ 1 $ 1 ) 2 + ^A 2 ($ 2 $ 2 ) 2 + A 3 ($i$i)($|$ 2 ) + A 4 ($I$2)(4 $ i) 

+ j^A 5 ($ t 1 $ 2 ) 2 + [A 6 (*t%) + A 7 (*J*2)](*}$2) + fc.cJ. (2.1) 

The coupling constants m^, m 22 and Ai 2,3,4 are real, while mi2 and As 6,7 can be complex 
if there are not any further restrictions. In the following we will however set them to be 
real such that there is no explicit CP-violation. 

The vacuum expectation values (VEVs) of $j are in general 

(*i)o 

and tan/3 is defined as the ratio of the tan/3 = v 2 /v\. 

The Higgs doublets can be rotated to a basis in which only one of the doublets has a 
vacuum expectation value using the angle /?. This is called the Higgs basis and is related 
to the general basis as 

Hi = cos @ $1 + sin ft e~ ie §2 , 

H 2 = - sin p $1 + cos P e~ ie ^ 2 , (2.3) 

















\ v 2 
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with 9 = 62 — 0\. Hence the VEVs for the doublets in the Higgs basis, with v 2 = v\ + v\, 
are 

c*>. = ^ 

We have defined both $j to have weak hypercharge +1. Doublets with weak hyper- 
charge — 1 can be constructed out of the complex conjugate fields via 

= i<72$* • (2.5) 

$1 and $2 consist of 8 real fields in total. Three of them correspond to the Goldstone 
bosons to be eaten by the weak gauge bosons and Z° upon spontaneous breaking of 
the gauge group SU(2)l x U(1)y- One of the standard conventions to write the doublets 
without the Goldstone bosons is (setting for clarity B\ = 0) 

-spH+ \ 
-^(cpv - s a h + c a H - is /3 A) J 

-^{spv + c a h + s a H + icpA) 




(2.6) 



Here H ± is the charged Higgs boson and the angle a (s a = sin a, c a = cos a) is introduced 
to diagonalize the CP eigenstates in the neutral sector, which can be divided into two CP 
even scalars: (H, h), and a CP odd pseudo-scalar: A. 

2.2 The Yukawa Sector 

The weak eigenstates of the SM fermions (with massless neutrinos for simplicity) are de- 
noted as 





The most general Yukawa interaction can then be written as 



(2.7) 



-C Y = Ql^Ur + Q^ivi D r + L L Q ir ftE R 

+Q L $2rttU R + Q l ^Dr + L L $ 2 rfiE R + h.c. . (2.8) 

We leave the generation index implicit here, all entities are matrices or vectors in the 
three-dimensional generation space. The rjf are the 3x3 matrices of Yukawa couplings 
for F = U,D, L. 
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In order to show more explicitly the physical content in the Yukawa couplings, we 
rotate the Yukawa coupling matrices to the Higgs basis by inverting Eq. (2.3) and inserting 
into Eq. (2.8). 

-C Y = Q L H lK u U R + Q L H lK gD R + LlH^Er 

+Q L H 2 pV U R + Q L H 2 p»D R + L L H 2P ^E R + h.c. . (2.9) 

The relations between the two sets of Yukawa matrices are 



4 

D 



cos Prft + sin/3(e-^), 
cos/37/f + sin/3(e +i0 77f ) , 

cos/3r/f + smf3(e +i9 r]%) ; (2.10) 



and 



$ = -smptf + cos/3(e-^) , 
p° = - sin /3r/f + cos /?(e +i % D ) , 

p% = -sin/3r/f + cos/3(e +ie r/ 2 L ). (2.11) 

The couplings to Hi produce the masses of the fermions. We can go over to the fermion 
mass basis by bi-diagonalizing the matrices k f with the unitary matrices V£, Vr'- 



K 



' - VEk$V? = ^Ml (2.12) 



v 



P F = V L F P^VP (2.13) 

The k f are diagonal, real and positive and are fully determined from the fermion masses 
M F with Mil = m u e t c - R F is still a general complex matrix whose non-diagonal matrix 
elements could cause tree level flavour-changing- neutral-currents. The reason is that we 
cannot in general diagonalize two different matrices simultaneously. The flavour changing 
charged currents are described by the matrix 

Vckm = VgV^ . (2.14) 

We now can derive the Yukawa interactions in the Higgs and fermion mass basis. 
Using the definitions of Eqs. (2.3), (2.6), and (2.10-2.13), the Yukawa interactions (2.9) 
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Type 


Ur 


Dr 


Lr 




P D 


P L 


I 


+ 


+ 


+ 


K U cot f3 


K D COt 


K L COt P 


II 


+ 






K U cot f3 


—k d tan P 


—k tan /3 


III/Y 


+ 




+ 


K U cot f3 


—k d tan /3 


k l cot P 


IV/X 


+ 


+ 




K U COt P 


k d cot P 


-k l tan/3 



Table 1. The different types of 2HDM with Z-i symmetry. The nomenclature follows [8]. The 
Z 2 charges for Higgs doublets are odd or —1 for $1 and even or +1 for $ 2 - The right-handed 
fermions have been given different Z2 charges assignment as shown. The Yukawa matrices p F are 
proportional to the k and thus also diagonal with the relation shown in the last three columns. 



become (see e.g. [14]) 



1 



K%_ Q + (^P R + p wl P L ) C/3 



Dt 



Dh 



V2 



K U S^ a + (p U PR + p U 'P L ) C p 



(p D P R + p w P L )sp. a \ DH + —D(p D P R - p w P L )DA 



Uh 



K U C!3- 



(p u P R + p u ^P L )s^ a ] UH - ±=U(p u P R - p u] P L )UA 



Lt 



K%_ a + {p^P R + P l 'Pl)C[} 



Lh 



K L C(3- 



(p L P R + p LJ P L )s^ a \ LH + ^L{P L PR ~ P L Pl)LA 



+ 



U(V C kmP D Pr ~ P U] V CKM P L )DH + + vp L P R LH + + h. 



c. 



(2.15) 



where Pr/l = (1 =t 7s)/2. One can clearly see, that if the Yukawa coupling matrices 
are not diagonal, there are flavour-changing-neutral-currents (FCNC) at tree level, which 
are absent in the Standard Model and are severely constrained by experiments. Therefore, 
either these terms are completely forbidden by certain symmetries or mechanisms, or they 
are sufficiently small to avoid the current experimental bounds. An early discussion is the 
paper by Glashow and Weinberg [2]. 

There are different known solutions to the FCNC problem. In this paper we study 
three different cases: 

• Z2 symmetry 

If there is only one Higgs doublet coupling to each type of fermions, the situation 
becomes the same as in the standard model. The FCNC couplings vanish completely, 
known as naturally vanishing FCNC [2]. An elegant way to achieve this is to impose 
a Z2 symmetry on the Lagrangian and set one of the Higgs doublets and some of the 
right handed fermions to be Z2 odd. The different cases depending on which fermions 
couple to the same doublets are listed in Table 1. We also note that the Higgs sector 
of the MSSM is of type II at tree-level. 
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Yukawa Alignment 

A more general way to diagonalize the Yukawa matrices simultaneously is the Yukawa 
Alignment model [3]. They proposed that the Yukawa coupling matrices rjf and r\ F 
are proportional to each other. So the rotated Yukawa coupling matrices k f and p F 
are also proportional to each other and can thus be diagonalized simultaneously. 

However, other than the models with Z2 symmetry, this alignment may be spoiled 
when going to different energy scales and some of the non-diagonal couplings leading 
to FCNC may become sizable. Studying limits on the proportionality constants from 
this source is one of the purposes of the present paper. 

Cheng-Sher Ansatz 

A third possibility is to keep the off-diagonal FCNC elements in the p F naturally 
small. The best known ansatz of this type was proposed by Cheng and Sher [9] 



P5 = AP^^. (2.16) 

The rrii are the different fermion masses. Since the diagonal elements of the k f 
have a hierarchy in size corresponding to the fermion mass hierarchy it is natural to 
introduce this also for the p F . The X F are expected to be of 0(1) and should be 
small enough to suppress FCNC to the observed level. We discuss these limits below. 
One should be aware that there are different parameterizations of the Cheng-Sher 
ansatz, some papers do not have the factor of \/2 in (2.16), e.g. [15]. 

2.3 RGE for Yukawa Couplings in 2HDM 

The variation of couplings and masses with the subtraction scale p is given by the renor- 
malization group equations (RGE). The running of Yukawa couplings in the 2HDM can be 
found in many places, e.g. [4, 10, 11]. We have also rederived them using the methods of 
[10]. 

Using the notation T> = 16-7T 2 d / '<i(ln p) the RGEs for the Yukawa couplings in the 
general basis are: 

2 

r ll: 



= -Aurt + £ Tr [N c («t + ^ ^t) + 

+\ E [vfv? + %V f ] f k + 4 E v?v? - 2 E [ 



= -An^ + £ Tr [N c (r£^ + fitf) + « 

+\ E + ] ^ + ^ E W - 2 E[' 



U U\ D 
% V k Vi 



-7- 



= -A m u k + E Tr [iV c ffirtf + i^i/f f ) + r^f 



E 



i 



L it i 1 L it L 



(2.17) 



where are given by the gauge couplings as follows 



A 



9 



17 



V 



(iV^-1) , f , x. , 
AT,. 53 + 4 52 + 12 51 



A D = - g[ 
A L 



15 2 9 2 
T 9i + ~ A 92 



(2.18) 



with gi = e/ cos0\y, 92 = e/ sin^^y, and g3 = g s , sinOw being the weak mixing angle. In 
turn the RGEs for the gauge couplings up to one loop level are 



V(g 2 ) 



i 



10 

+ 7T n ? J 9i 



n q 9 2 



■~{llN c -2n q )gl. 



(2.19) 



n q is the number of active quarks above energy threshold. In this paper we will always use 
n q = 6. We have checked that using the two- loop running for g% produces only a small 
change in our results. We thus expect that the effect of running with n q = 6 from mz to 
the top threshold rather than n q = 5 will not introduce a significant effect. 

Finally the RGEs for the fields and thus for the vacuum expectation values e *Vi are: 



V(e i9k v k 



+ [\al + \ 9 l ) e^v k 



e We v £ 



(2.20) 



Note that the running of the Yukawa couplings as given in (2.17) is independent of 
the couplings in the Higgs potential (2.1). They only appear at the two-loop level. Thus 
we limit ourselves to studying the evolution of the Yukawa sector by itself and do not 
include the evolution of the parameters of the Higgs potential. One should keep in mind 
that the evolution of the latter could also signal the breakdown of a given model. This 
has for example been studied in [16] although only including the top Yukawa coupling. A 
complete one-loop treatment of the Higgs sector would require also the inclusion of the 
complete Yukawa sector at one- loop. We foresee to include this in future versions of the 
2HDMC calculator [17]. 

Using the definitions (2.3), (2.10) and (2.11), the RGEs can be rewritten in the Higgs 
basis. The vacuum expectation value v, the phase difference between the two vacuum 
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expectation values 6 and the angle j3 relating the general basis and the Higgs basis satisfy 
the following RGEs: 



V {v 2 ) = -2Tr [iV c + k£ + 44* 



V 2 + 



3 2 9 2 ' 
2»i + 2* 



£>(tan/3) 



2 cos 2 



Tr 



^ (pW + «W f + "fob* + /^o Dt ) 



+ K oPo^ + Po ^O^ 



isin(2/3) 



Tr 



( L L\ L L}\ 

y K o Po Po K o j 



(2.21) 



Finally the Yukawa couplings in the Higgs basis, in other words the matrices Kq and Pq 
satisfy: 



V {4) = -A UK % + Tr [iV c [4 4 + + 4^0 



-1 tan/3Tr(iV c _ ^t) _ Ng ( K ^t _ 



4 f ) 



+ < 



PoVo * + PoVo ^ + K K ^ + ^0*^0 ^ 



ft U i t/f [/ 

Po Po + K o K o 



-2PqKq ^ Pq — 2KqKq ^ Kq > , 



(2.22) 



Kq 



V (4) = -A DK g + Tr [N c (4 4 + ngn^) + 

+ \ tan/3Tr(iV c (^t _ ^t) _ Nc ( K ^t _ p D ^ 



( L L] L L\\ 1 D 

y K oPo Po K o J ( K 



+ < 



Po^Po * + PoVo * + K o'' c o * + K o K o ^ 



K o* + K o* 



D\ D , £»t £> 
Po Po + K K 



-2pqK q ^ Pq — 2kqK q ^Kq > , 



(2.23) 



+ ^tan/3Tr(iV c _ ^t) _ Ng (tftf _ p D ^ 



L L\ L L\\ I L 

K o Po ~ Po K o I ( K o 



+ 



1 ( L L\ , L L]\ L . L ( L\ L , \L l\ 

2 \^Po Po + K o K o J K o + K o ^Po Po ~i~ K o K o J > 



(2.24) 



^ + 2Tr iV c ( ^ ^ + | + ^ 



Kg 



+Tr 
1 



Po 7 



cot/3Tr|iV c _ ^t) _ Nc ^dd, _ p o K D^ 



( L L\ L L\\ 

y K oPo Po K o J 



>P U o 



+< 



PoPo ^ + Po Po ^ + K o K o ^ + K K ^ 



P U o+P U o 



Po Po + K o K o 



~2po Po ^ Po — 2^0 Po ^ K o 



(2.25) 



V(p°) = -A DP o° + 2Tr [iV c (k% + + p^ 

\n,{p U opT + P^P D o ] + P L oP L o ] ) 



K o 



+TY 
+ 



P°o 



~~ ( K o Po^ ~~ Po K > 



f f ) Po D 



+< 



2 L 



Po Po + Po Po K K K K 



Po^ + P? 



Po Po ~r K K 



(2.26) 



- 10 - 



Vpl = -A L p% + 2Tr|iV c (4p U ^ + p° k^) + p^ f J>4 



■ N, 



y K o Po 



Po D -o Dt ) 



L L\ 
K Po 



L L\ 
Po K 




+ -. 



L L\ . L L\ 
Po Po ' K K 



Po + Po (poVo 



4- k AL k L 



(2.27) 



Before ending this section we note that the tan f3 dependent terms in the evolution 
equations for the Yukawa couplings disappear in the real case. In the CP-violating case p 
is no longer basis-independent and therefore there is a residual dependence on tan (3 in this 
case. For a thorough discussion of basis independent quantities in the CP-violating case 
we refer to [18]. 

3 Constraints and SM input 



3.1 



Low-energy constraints on \f- 



In the recent review of 2HDM [11], the authors have given a comprehensive overview on 
the latest constraints on the \fy The most stringent ones are in the quark sector, coming 
from the neutral meson mixing, and we will therefore limit ourselves to these constraints 
in the following. 

The master formula for F° — F° mixing mediated by tree level Higgs scalars in the 
vacuum insertion approximation can be found in [19]: 



AM F 
S ^ 
Pf 



\P F \ 2 
M F 



S F 



B F f F M 2 F 



1 

6 

\B F f F M 2 F 



mr 



1 + 



+ 



+ 



Pf 



H 

Ml 



m , 



(3.1) 



1 + 



(mi + rrij) 
llMf 



(mi + m 



j > 



Here Mp and AMp are the mass and mass difference of the neutral mesons respectively, 
and jp is the corresponding pseudo-scalar decay constant. The parameter Bp is defined 
as the ratio of the actual matrix element compared to its value in the vacuum insertion 
approximation [19]. The numerical values of the parameters we use are listed in Table 2. 

To calculate the limits on , we require that the sum of the SM and 2HDM theoretical 
predictions for AMp does not exceed the experimental value by more than 2 standard 
deviations: 



AM F M + AM| HDM < AM° xpt + 2a 



(3.2) 
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Meson 


M F (GeV) 


B F 


If (GeV) 


K° (ds) 


0.4976 [20] 


0.75 ±0.026 [21] 


0.1558 ±0.0017 [21] 


D° (uc) 


1.8648 [20] 


0.82 ±0.01 [22] 


0.165 [22] 


B° d (db) 


5.2795 [20] 


1.26 ±0.11 [21] 


0.1928 ±0.0099 [21] 


B° s (sb) 


5.3663 [20] 


1.33 ±0.06 [21] 


0.2388 ± 0.0095 [21] 



Table 2. Parameters of the neutral mesons K°, D a ,B^ and B a s . 



where a = y o^xpt + °sm ^ s a combination of the experimental and theoretical uncertain- 
ties. For the K° - K° and D° - D° mixing, the non-perturbative interactions make the 
SM calculation very difficult. Here we therefore simply assume that the 2HDM contribu- 
tion is not larger than the experimental value by more than 2 standard deviations. This 
corresponds to setting the SM contribution to zero in Eq. (3.2) as was done in [15]. The 
experimental and SM values we thus use are listed below. 



1. K° - K°: 



AM^ x pt = (3.483 ± 0.006) x 10~ 15 GeV 
AM|M = 



2. D° - D°: 



3- B° " B°: 



AM^f = 1.57j#4i5 x 10~ 14 GeV [20] 
AMfM = 



AM£ pt = (3.344 ± 0.0197 ± 0.0197) x 10~ 13 GeV [20] 
AM| M = 3.653±{j;io x 10- 13 GeV [23] 



4. B° - B°: 



AM^ pt = (116.668 ± 0.270 ± 0.171) x 10~ 13 GeV [24] 
AM|f = HO.6^^9 1 x 10~ 13 GeV [23] 

The 2HDM contribution is then calculated using Eq. (3.1). We note that the quark 
masses appearing in Eq. (3.1) are the low energy ones defined more or less at the scale of 
the respective meson masses. For internal consistency we use the following values from ref. 
[25] (in GeV): 

m u (2 GeV) = 2.2 x 10~ 3 , m c (m c ) = 1.25 ; 

m d {2 GeV) = 5.0 x 10~ 3 , m s {2 GeV) = 0.095 , m b (m b ) = 4.2 . 
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However, the impact of the actual quark masses used is very small since the masses ap- 
pearing in \pfj\ 2 and the dominant pseudo-scalar matrix element Mp/(mi + rrij) 2 essentially 
cancel, and we get similar results using the masses defined at mz instead. 

From Eq. (3.1) we can see that the main uncertainty of this estimate is due to the 
unknown masses of the CP-even and CP-odd Higgs bosons. It is also clear that the con- 
tribution to the mixing from the CP-odd exchange is much larger due to the extra factor 
11 in Pp for the dominant pseudo-scalar matrix element. We will consider three different 
representative cases. We also remind the reader that in some cases there is no factor of \pl 
in the definition of Xfj. With all this in mind we get the following constraints on A^: 

• rrih = tuh = rriA = 120 GeV 

Kc < 0.13, 

X ds < 0.08, X db < 0.03, X sb < 0.05 . 

• rrih = run = mA = 400 GeV 

Xuc < 0.44, 

X ds < 0.27, A db <0.12, A s6 <0.18. 

• m h = m H = 120 GeV m A = 400 GeV 

X uc < 0.30 , 

A ds <0.20, A df) <0.08, A sfe <0.12. 

The first and second cases are examples of typical low and intermediate masses for 
the Higgs bosons, whereas the last case illustrates that the main restriction comes from 
the exchange of the CP-odd Higgs. All in all we conclude from these different cases that 
a representative value for these constraints is given by A£. - < 0.1 and this is the generic 
value we will use when analyzing the effects of Z2 breaking in the running of the Yukawa 
couplings in the next section. 

3.2 General input 

For the RGE evolution towards high scales we need a set of input parameters at the low 
scale fj, = mz = 91.186 GeV. The experimental input we have are the masses and the 
measured parameters of the CKM-mixing matrix as well as the gauge couplings. We have 
neglected constraints coming from the neutrino sector. The quark and charged lepton 
masses at the scale mz we take from Ref. [25], their values are (in GeV) 

m u = 1.29 x 10~ 3 , m c = 0.619, m t = 171.7; 
m d = 2.93 x 10 -3 , m s = 0.055, m b = 2.89 ; 
m e = 0.487 x 10~ 3 , m„ = 0.103, m T = 1.746. 
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For the 3x3 CKM matrix we use the PDG [20] phase convention 



Vckm 



C12C13 S12C13 size~ lS 

-S12C23 - ci2S 2 ?,size l5 C12C23 - si 2 S2?,size l5 S23C13 I • (3.3) 

\ S12S23 - Ci 2 C23S ri e lS -C12S23 - S 12 C23Sl3e lS C23C13 



where Sij = sin 9ij and Cij = cos 9ij . We will also use this convention for the phases at the 
high scale. The values for the angles and the phase follow from [20] 

S21 = A , S23 = AX 2 , 
iS _ A\ 3 (p + iy)Vl-A2\* 
13 v / T^A2[l-^ 2 A 4 (p + if/)] ' 

with 

A = 0.2253, A = 0.808, p = 0.132, fj = 0.341 . (3.5) 

There is of course still a large freedom in how one chooses the remaining freedom at 
the weak scale mz- We chose to put the CKM-mixing always in the down quark sector 
and have thus at the EW scale 

vE = vK = i 
vt = vk = i. 

The last two are a consequence of our neglecting neutrino masses and mixings. The Yukawa 
couplings at the EW scale are thus: 

( K o)ij = K fj = i (Po)ij = Pij = U,C,t) 



(«o )ij = Vqkm = Vckm , (p )v = ^ckm ftj = d, s, b) 



At any energy higher than the EW scale, the Yukawa couplings kq and po in general 
become non-diagonal and complex. Thus they need to be transformed to the mass eigen- 
states by the bi-diagonalization defined in Eq. (2.13) in order to give k and p. The latter 
can then be used together with the diagonal elements of the former to calculate • . When 
performing the bi-diagonalization we always keep to the PDG conventions for how to write 
the CKM matrix. 

For the electroweak VEV we use v 2 = l/(y/2G F ) with G F = 1.16637 • 10~ 5 GeV" 2 
from PDG [20] and for the phase difference between the two VEVs we start from 9 = such 
that there is no spontaneous CP-violation. For the gauge couplings we use the PDG [20] 
values: a = 1/127.91, a s = 0.118 and for the weak mixing angle we use the on-shell value 
sin 2 9 W = 0.2233. 
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Type 



II 

III/Y 

IV/X 



I 



l/tan/3 
l/tan/3 
l/tan/3 
l/tan/3 



l/tan/3 

— tan /3 

— tan /3 
l/tan/3 



1 / tan /3 

— tan /3 
1 / tan /3 

— tan (3 



Table 3. The diagonal A£ 



in 2HDM models with Z 2 symmetry. 



4 RGE analysis 

We have implemented the RGE equations in the Higgs basis given above in three different 
computer codes. The matrix operations have been performed with either the C++ template 
library Eigen [26] or the GNU Scientific Library (GSL)[27] and the in total 114 ordinary 
differential equations are handled by the ODE-solver in GSL using the explicit Runge- 
Kutta-Fehlberg (4,5) method. The programs have been tested against each other and also 
by comparing with the results from [10]. 

In this section we will start by briefly exploring the behavior of ^-symmetric models 
and then study a number of ^-breaking models in more detail. 

4.1 .^-symmetric models 

From table 1 and the definitions of k f and p F , we get the diagonal elements of A^ in 
terms of tan/3 for the four different 2HDM types as shown in table 3. Since in this case 
the Yukawa couplings are given by tan /3 it is a real physical parameter. In addition the 
evolution of the Yukawa couplings will only depend on the initial value of tan /3. 

Since the ^-symmetry is enforced the Yukawa couplings stay diagonal and the only 
thing that can happen during the evolution is that one or more of the Yukawas will blow 
up due to the presence of a Landau pole. This signals the breakdown of the perturbative 
description and calls for a new theory at the corresponding energy scale. The position of 
the Landau pole will depend on the initial value of tan /3 and which of the four types we 
are considering. 

In Fig. 1 we show the position of the Landau pole as a function of the input tan/3. 
For the lower limits, the results are almost the same for all types, and the lines are more 
or less on top of each other. This is natural since in this regime the evolution is essentially 
driven by Xu, which is the same in all types. For the upper limits, on the other hand, 
there are some differences. First of all there is no upper limit on tan /3 in the type I model, 
which means there is no Landau pole below 10 16 GeV if the input tan/3 > 1.1. For the 
other types, the upper limits are shown in the right panel of Fig. 1. The differences can be 
understood from whether the evolution is driven by A^ (type III), A TT (type IV) or both 
(type II). 

4.2 .^-breaking models 

Before starting to analyze the ^-breaking models we note that, as shown by [4], the Z2- 
symmetry of the RGE's is still preserved if all the X F, s for the different types are rescaled 
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Lower limits of tanfS 



Upper limits of tanfi 




Figure 1. The starting value of tan/3 as a function of the position of the corresponding Landau 
pole (A 

Landau— pole) hi the different 2HDM types. There are lower limits on tan/3 for all four types 
(left), but only type II, type III/X and type IV/Y have an upper limit of tan/3 (right). 



with a factor x for the cot (5 ones and 1/x for the tan/3 ones. In addition, when tan/3 is no 
longer related to the Yukawa couplings it does not have any physical meaning, since it only 
reflects the basis choice for the general 2HDM. In the following we will only be considering 
cases with p real at the starting scale. This means that the only source of CP-violation 
is from the CKM-matrix. Thus the CP-violating effects will be small and therefore the 
dependence on tan /3 very limited. We have verified this numerically for a number of cases 
and in the following we set tan/3 = 1. 

In this subsection, we will also explore the non-diagonal elements of FCNC Yukawa 
couplings. We know that in the Z2 symmetric case, the tree level FCNC couplings will 
remain equal to zero (up to the numerical precision) up to arbitrarily high energy scales 
since they are protected by the symmetry. However once we break the Z2 symmetry in 
some way, this protection is not effective anymore and the off-diagonal elements Xf,- may 
start to grow. 

The actual values of the non-diagonal FCNC Yukawa couplings pf\j at different energy 
scales will depend on how much we break the Z2 symmetry. We can thus use the size of the 
X^j.j as a measure of how severe different types of Z2 symmetry breaking are. Of course 
we do not know how large the AEy can be at higher scales. Still it is reasonable to assume 
that the values should not be widely different from the ones at the EW scale. Thus we will 
use a generic value of Ajy < 0.1 as a limit on how much Z2 symmetry breaking should be 
allowed and see at which energy scale this limit is reached. 

The argument behind this is essentially that we can use the RGE evolution to analyze 
the stability of the assumptions underlying different 2HDMs under variations of the scale 
where the model is defined. A large sensitivity indicates that the assumptions behind the 
model are not stable meaning that they are either fine-tuned or incomplete such that there 
for example will be additional particles appearing when going to a higher energy. From 
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this respect we will thus study both the appearance of a Landau pole as well as off-diagonal 
Yukawa couplings leading to FCNC larger than experimentally allowed at the EW scale. 
We also note that as will become clear below there is a small dependence on at which 
scale we apply the above argument. Requiring stability up to 10 3 GeV gives very similar 
constraints on the amount of ^-breaking that is allowed as when using 10 15 GeV. 

As an alternative way of assessing the amount of .^-breaking that is allowed by the 
experimental constraints from FCNC at the EW scale we have also considered a set-up 
where the Z2-symmetry is broken at a high scale. In this set-up we start from a Z%- 
symmetric model at the EW scale and then evolve it to the high scale of interest. The 
resulting model is then used as the starting point for exploring different ways of breaking 
the Z<i symmetry. Once the .^-breaking has been introduced the different models are then 
evolved down to the EW scale for comparison with the experimental constraints. We have 
verified that in representative cases the results obtained in this way are very similar to the 
first approach and therefore we will not go into any more details. 

There are many possibilities to break the Z2 symmetry and in the following we will 
consider three ways: aligned, diagonal and non-diagonal Xfj as defined below. In most 
cases we will concentrate on the effects of breaking the symmetry starting from a type I or 
type II model. The reasons for this is on the one hand that these models are the most well 
studied cases in the literature and on the other hand that it is in the quark sector that we 
have the most stringent constraints on the FCNC Yukawa couplings. Thus the breaking of 
the Z<i symmetry in the lepton sector will typically have small effects. 

In order to be able to separate the effects of breaking the Z<i symmetry in different 
ways we will limits ourselves to breaking the symmetry in one specific way at a time. 

We start by noting that in the Z2 symmetric models at least two of the A^ are always 
equal whereas the third one is the same as the other two in type I and the negative inverse 
of them in the other types. When going to the aligned models we will therefore keep two of 
the sectors in fulfillment with the ^-symmetry and only break the symmetry through the 
relation to the third sector. In other words either setting X® = A^, A^ = Xf i: or A^ = A§ 
and letting Xf { of the third sector vary independently of the other two. 

Another way of breaking the ^-symmetry is by keeping the A^ diagonal but letting 
the individual diagonal elements be non-equal as has been studied by Mahmoudi and 
Stal [8]. We will analyze the effects of this type of breaking in the up and down sectors 
separately again starting from the ^-symmetric cases with either X 1 ? = A^ = Xu or 
A|? = A^j = — 1/Xtt- In other words using the type I or II ^-symmetries as starting point. 

The third way of breaking the ^-symmetry that we will consider is by setting the 
non-diagonal elements of A^ nonzero already at the starting scale. Again we will consider 
setting the up-sector and down-sector non-diagonal elements non-zero separately and apply 
the type I or type II symmetries for the diagonal elements. 

4.2.1 Aligned models 

We start by analyzing the three different versions of Aligned models with A^, X D , and X L 
pairwise equal. Based on the similarities with the /^-symmetric models we call them I /II, 
III, and IV respectively and their free parameters are as follows 
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Landau pole ^. D = X L 



Landau pole X u = X L 
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Figure 2. The energy scale at which the Landau pole is encountered as a function of pairwise 
combinations of the starting values for A^{, and \f { as indicated in the figure for the three 
different versions of aligned models explained in the text. The areas inside a given contour are 
allowed by the requirement of not having a Landau pole. The different contours are as follows 
starting from the center: 10 15 , 10 10 , 10 5 , 10 3 , and 300 GeV. 



. Aligned I/IL X%, \° = A|j 
. Aligned III: \g, \V. = \L 
. Aligned IV: A|, Ag = \? 

First we consider the effects of requiring that there is no Landau pole encountered 
when evolving to higher scales. We therefore plot in Fig. 2 the scale at which the Landau 
pole is reached as a function of the starting values for pairs of X u , X D , and X L . This means 
that for a given energy scale the points inside the corresponding contour is allowed by this 
requirement. As can be seen from the figure, the position of the Landau poles is very 



- 18 - 



X D = X*~ or %y= X 





i 




: ;> 5 1 1 y 


4 






/ / /» i; \ ' /_ 

■ |: '• jf 

//Ik v /f 


2 






• / -y / >'•. """ / :' 300 - 

/.■'/ X ■ 






. \ 
\ \ 




-2 




\ 

■■ / j 




-4 


/ i 




\|/t/ / / 
ills / i i 



-4 -2 2 4 



Figure 3. Same as Fig. 2 but also applying the constraints from the non-diagonal Xf^.j. The plot 
shows the results for = A^, but the same results are also obtained for Afj = A^. 

similar to the situation for the ^-symmetric cases and there is only a small correlation 
between the values of the aligned Yukawas where the Landau pole is reached. 

Applying also the condition that the off-diagonal elements of should respect the limits 
given by the meson mixing constraints also at higher scales has a potentially large impact 
on the allowed regions. This is the case for the aligned models of type I/II and III, where 
X L is set equal to X D and A^ respectively, as can be seen in Fig. 3. In fact, within the 
parameter region displayed in the figure (note the difference in scale compared to Fig. 2) 
there is no difference between the two cases and therefore we only show one of them. 
However, as may also have been expected, there are no additional constraints in the case 
when A^ and A^ are set equal since the off-diagonal lepton Yukawas are always small as a 
consequence of the small lepton masses and the limited cross-talk between the quarks and 
leptons. In other words breaking the Z2 symmetry between the quarks and leptons has no 
effect in this respect. 

For reference we have also included lines corresponding to the Z2 symmetric relations 
in Fig. 3. Along these lines it is the Landau pole that gives the limit but in the other regions 
the limit comes from the off-diagonal elements. We also note that the plot is symmetric 
under inversion through the origin (x,y) — > (—x,—y), which follows since the evolution 
equations for are all even under — > —p^ as long as the imaginary parts of Kq and 
Pq are small. 

It is also interesting to compare the results for non-equal X D and X u with the con- 
straints on Xbb and Xu obtained from b — > 57 in [8]. Applying the conditions of stability 
when evolving to higher scales and that the non-diagonal Yukawas should stay small es- 
sentially removes the regions |A^| > 1 including the fine-tuned regions where Xbb and Xu 
are both large (> 2) and have the same sign. 
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Figure 4. The constraints on the starting values of £ = A^/A^ (left) and £ = — A^A^ (right) as 
a function of the renormalization scale where the off-diagonal elements reaches 0.1 in the Aligned 
models of type I and type II respectively for the representative values Xf t = 0.02 and 0.5. 



As special cases we also show in Fig. 4 the results for A^ = 0.02, 0.5 and either X® = 
Xf { = £A^ (type I) or A§ = Xf { = — £/A^ (type II). From these plots it is clear that for 
A^ = 0.5, the off-diagonal elements puts strong constraints on the ^-symmetry breaking 
parameter £ = X^/X^ ( £ = — X^Xf { ) with typical values being £ < 3 — 10 (2 — 5) for type 
I (II). For A^ = 0.02 on the other hand the constraints are very mild in a type I set-up 
with £ < 100 — 1000 allowed, whereas in a type II setup only £ values very close to 1 or 
£ < 0.05 — 0.1 are allowed. The two possibilities corresponds to two distinct regions in the 
A^,A^ plane. The first one where X® « — 1/A^ and the second one where X® is small 
(< 2 — 5). For comparison we recall that the Landau pole constrains A^ < 70 — 200 more 
or less irrespectively of A^. So the constraints on £ are more or less trivial in this case. 

4.2.2 Diagonal models 

Next we consider in more detail models with .^-breaking in either the up or the down 
sector. To make the discussion more clear we only consider models where A# and A^ are 
related in a Z2 symmetric way and since we have seen that the effects of the lepton sector 
is small we always set A^ = A^. (If Xu and are not related in a Z2 symmetric way then 
we are more or less back in the aligned models since these two are the dominant Yukawas). 
In other words we only partially break the alignment. 

Thus we start with considering ./^-breaking in the up-sector with X D = Xu (type I) or 
X D = —1/Xtt (type II). For simplicity we also set X uu = A cc . 

First of all, as we show in Fig. 5, the Landau pole gives the restriction A cc < 400 — 500 
both for type I and II, again more or less independently of the value of Xu- We also want 
to emphasize that even though it is not really discernable from the figure, there is also a 
lower limit on Xu > 0.01 from the Landau pole for Ay, for type II. 

The figure also shows that the impact of constraining the off-diagonal elements to be 
less than 0.1 is limited for the type I set-up. In fact for Xu = there is not additional 
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Figure 5. The energy scale where the Landau pole is reached (upper panels) together with the 
scale where one of the non-diagonal A|^ ■ = 0.1 (lower panels) as a function of the input values A cc 
and At*. In the left (right) panels X D = X L = Att(— 1/Att). 



constraint from the off-diagonal elements. In the type II set-up the constraints are more 
severe but even so quite mild. 

To get a better picture of the range of the amount of ^-breaking allowed we also give 
in Fig. 6 the constraints on the ratio A cc / Att in type I and II set ups for our standard values 
Att = 0.02 and 0.5. From the plots it is clear that this ratio can be as large as ~ 1000 
without generating off-diagonal A^ > 0.1 all the way up to the GUT scale. 

Next we consider ^-breaking in the down-sector with A^ = A^ (type I) or Ay, = 
— 1/ A Yi (type II) . Similarly to the up-sector we set \m = X ss for simplicity. Also in this case 
the constraints from the Landau pole are similar for the two set-ups with \ ss < 400 — 700 
in both cases with a small correlation with the value of A^ and A^ for a set up of type 
I and type II respectively as can be seen from Fig. 7 (upper panels). However, contrary 
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Figure 6. The constraints on the input values £ = X cc /Xu as a function of the renormalization 
scale where the off-diagonal elements reaches 0.1 in the diagonal models of type I (left) and type II 
(right) for the representative values Xu = 0.02 and 0.5. 



to the up-sector the figure (lower panels) also shows that the effects from requiring the 
off-diagonal Yukawas to be small are quite severe. In the type II case one can even see a 
mild preference for solutions with X ss ~ A55. 

To get a more quantitative picture of the constraints we show in Fig. 8 the ratio 
£ = Xss/Xbb for type I and type II using the values A^ = 0.02 and 0.5. In the type II set-up 
the constraints are especially restrictive with £ < 4 — 10 for A-{ = 0.02. In the type I set-up 
the constraints are less severe but even so stronger than the corresponding ones from the 
up-sector. 

4.2.3 Non-diagonal models 

Finally we consider the case of breaking the ^-symmetry from having non-zero non- 
diagonal elements in the up- or down sectors. As starting point we again use the Z2 
symmetric models of type I or II for the diagonal elements and then set either = 0.1 
or ASy = 0.1 at the EW scale in order to break the Z2 symmetry. 

Quite unexpectedly the additional constraints from requiring the off-diagonal elements 
to stay small are limited. The corresponding plots for the case of only considering the 
Landau pole are essentially straight vertical lines. Thus we do not show the effects of 
applying the two constraints separately. In fact it is only in case II with A|? = Ajj = — 1/A^ 
and Xf,j = 0.1 that the requirement of having Xf_^j(fi) < 0.1 gives any discernable effect 
and then only for small X^ < 0.2. On the other hand, in this case the constraints are 
very strong as also illustrated in Fig. 10. It is interesting to note that it is actually the 
off-diagonal elements in the down-sector that become large whereas the ones in the up- 
sector remain in accord with the limit A^ -(//) < 0.1. This means that even though there 
are presently no direct experimental constraints on X c t and X u t they are in this case highly 
constrained from the link to the down-sector through the RGE evolution. This is then the 
case in the MSSM, the prime example of a type II 2HDM, for large tan /3. To see more 
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Figure 7. The energy scale where the Landau pole is reached (upper panels) together with the 
scale where one of the non-diagonal X 1 ?,. =0.1 (lower panels) as a function of A ss and A^. In the 
left (right) panels = \bb{— 1/Xbb) and in all cases X L = Xbb- 



clearly what happens we show also in Fig. 10 the RGE evolution of the relevant off-diagonal 
elements for the input values A^ = 0.02 and XV,- = 0.001, A£. ■ = 0. 

5 Conclusion 

We have seen that the RGE evolution is a useful tool to analyze the stability of the as- 
sumptions underlying different versions of the 2HDM under variations of the scale where 
the model is defined. A large sensitivity indicates that the assumptions behind the model 
are not stable meaning that they are either fine-tuned or incomplete such that there for 
example will be additional particles appearing when going to a higher energy. From this re- 
spect we have studied both the appearance of a Landau pole as well as off-diagonal Yukawa 
couplings leading to FCNC larger than experimentally allowed at the EW scale. 
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Figure 8. The constraints on the input values £ = X ss /Xbb as a function of the renormalization 
scale where the off-diagonal elements reaches 0.1 in the diagonal models of type I (left) and type II 
(right) for the representative values Xu = 0.02 and 0.5. 



Based on our studies we have seen that the constraints from avoiding a Landau-pole 
are in general the same irrespective of the ^-symmetry. They appear as soon as the 
magnitude of one of the Yukawa couplings becomes of order 1. 

The constraints from the off-diagonal elements on the other hand depend on the details 
of how the ^-symmetry is broken: 

• breaking the Zi relation between X D and X u as in the Aligned models is highly 
constrained with X D /X U < 10 or -X D X U < 10 unless X D and X u are both < 2, 

• breaking it instead in the up-sector by having A cc and Xu non-equal gives a small 
difference compared to the constraints coming from the Landau pole with ratios 
A cc /Att > 1000 allowed, 

• in the down sector the constraints can be much stronger, but also more dependent 
on the relation between Ay, and Xu, ranging from X ss /X bb < 10 for X bb = 50 and 
X tt = -0.02 to X ss /X bb < 10000 for X bb = X tt = 0.02, 

• in the non-diagonal case the constraints are weak when starting from A^- = 0.1 
except in the case X bb = —50 and Xu = 0.02 where we find AEy < 0.001. In all cases 
it is the Xf,- that become large. 

From this we can conclude that starting from a type I Z2 symmetry there is quite a 
lot of room for breaking the symmetry as long as one does not encounter a Landau pole 
except that X D /X u < 10 has to be respected. In the type II case however, the room for 
breaking the symmetry is much smaller for large X bb . This is natural since in the latter 
case Ku and p bb are both large. It is also interesting to note that this corresponds to the 
situation in the MSSM with large tan (3. Finally we conclude that there is little hope to 
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Figure 9. The constraints from the Landau pole and the off-diagonal elements as a function of A^ 
and the off-diagonal elements A^- (up) or A^y (down) at the input scale for the type I (left) and 
type II (right) relations for the diagonal elements. 



see effects of non-diagonal Yukawa couplings in the top system in a type II model such as 
MSSM if tan/3 is large. 
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Figure 10. Left: The starting value Xf,, as a function of the energy scale \x where one of the non- 
diagonal elements of \V£ becomes larger than 0.1 for Xf i — 0.02 and a type II relation between the 
diagonal elements, A® — \\ { — — 1/A^. Right: The RGE-evolution of the non-diagonal elements 
A c t(/x), A s f,(/i), and \db(n) in the same case for AZy = 0.001. 
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